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In this paper we analyze the possibility of persistent currents of a two-species bosonic mixture 
in the one-dimensional ring geometry. We extend the arguments used by Bloch [lj to obtain a 
criterion for the stability of persistent currents for the two-species system. If the mass ratio of the 
two species is a rational number, persistent currents can be stable at multiples of a certain total 
angular momenta. We show that the Bloch criterion can also be viewed as a Landau criterion 
involving the elementary excitations of the system. Our analysis reveals that persistent currents at 
higher angular momenta are more stable for the two-species system than previously thought. 
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I. INTRODUCTION 

The hallmark of superfluidity is the possibility of dissipationless flow in situations where the flow of a normal fluid 
t>X) would degrade as a result of viscosity. The textbook example of this is the flow of a superfluid through a narrow 
capillary [2| ■ According to the Landau criterion [3| , the superfluid component flows without dissipation provided the 
superfluid velocity does not exceed some critical value. In this situation, the normal component remains locked to the 
walls of the capillary whereas the superfluid component, carrying zero entropy, flows as if the walls of the capillary 
behaved as a perfectly smooth conduit. If the capillary is now bent into a torus, one can imagine that a flow, once 
established, could persist indefinitely. 

The conditions under which persistent currents can occur for a bosonic mixture in the ring geometry is the subject of 
this paper. The usual analysis Q leading to the Landau criterion is not obviously applicable since one cannot invoke 
Galilean invariance for this closed system. However, for a system having a single component, Bloch |l| presented 
general arguments based on an analysis of the quantum mechanical many-body wave function which provided a 
criterion for persistent currents. He considered an idealized one-dimensional ring of radius R in which the particles 
interact via an arbitrary pair-wise interaction. Since the total angular momentum commutes with the Hamiltonian of 
the system, the stationary states have energies E a (L) which are functions of the angular momentum quantum number 
L; all other quantum numbers are subsumed in the index a. Bloch showed that these energy eigenvalues take the 
form 

>: 

where Mr = NM is the total mass of the system containing N particles of mass M. The first term on the right hand 
t"*** ' side of Eq. ([TJ is interpreted as the kinetic energy of a rigid ring rotating with angular velocity f2 = L/MtR 2 ■ The 
second term, e a (L), corresponds to internal excitations of the system; it has the periodicity property 



e a (L + Nh) = e a {L). (2) 

This implies that the system can find itself in the same internal state for angular momenta that differ from each other 
by multiples of Nh. In addition, e a (L) has the inversion property 

e a (—L) = e a (L), (3) 

which reflects the fact that the energy does not depend on the sense of the angular momentum. 

The state with the lowest energy for a given L will be given the label a = 0. In the noninteracting limit, eo(L) 
has a local minimum at L = [l|; one expects this property to persist with repulsive interactions. The periodicity of 
this function then implies that Eq(L) can exhibit local minima at certain multiples of Nh. If and when such minima 
occur, Bloch argued that the system is capable of sustaining persistent currents. Conversely, if Eq(L) is not at a 
local minimum, nonidealities will induce transitions which change the angular momentum and hence the flow of the 
superfluid around the ring. 

In Sec. HIl we extend Bloch's analysis to a two-species gas containing Na particles of type A and Nb particles 
of type B. Here the term "species" can refer either to different kinds of atoms or to atoms distinguished by their 
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hyperfine states. When the masses of the two species are different, we find that the energy can still be written in the 
form of Eq. ([!} but in general, eo(L) is no longer a periodic function of L. However, if the masses are equal, 6q(L) is 
found to have the same periodicity as for the single-species case with ./V = Na + Nb ■ In the case that the mass ratio 
Ma/Mb is a rational number, eo{L) remains a periodic function of L but with a periodicity that differs from Nh. For 
these special cases, Bloch's arguments for the possibility of persistent currents goes through as for the single-species 
case. For arbitrary mass ratios, Eq(L) may still exhibit a local minimum at sonic finite value of L but there is no 
general argument which can be used to determine where such a local minimum might occur. 

We go on to show that Bloch's criterion for persistent currents can be phrased in terms of the more familiar Landau 
criterion. For Ma = Mb, eo(L) is periodic and a Landau criterion can be formulated at the discrete set of angular 
momenta L = L n = nNh, with n an integer, where the system can be taken to be in its internal ground state. The 
Landau criterion then imposes a constraint on the spectrum of the elementary excitations with angular momentum 
rati and energy e(m). In Sec. IIII1 these excitation energies are determined for the two-species system in the Bogoliubov 
approximation. In general there are two Bogoliubov modes which are usually phonon-likc at long wavelengths. For 
the case Ma = Mb, the Landau criterion then suggests that persistent currents may be stable for certain values of n. 
However, if the interaction parameters satisfy a certain relation (given in Sec. IHIj) . one of the Bogoliubov modes has 
a particle-like dispersion and the Landau criterion leads to the conclusion that persistent currents are unstable for all 
n. 

The above conclusion was arrived at earlier by Smyrnakis et al. Q based on an analysis of the mean-field Gross- 
Pitaevskii (GP) energy functional for the two-species system. With the assumption that all interaction parameters 
are equal, these authors determine Eo(L) by minimizing the GP energy functional subject to the constraint that the 
average angular momentum of the system is L. Although persistent currents are destabilized at L = L n , the authors 
find that Eq(L) can exhibit local minima at non- integral values of I — L/Nh. In particular, they show that persistent 
currents are stable at I = xa = Na/(Na + Nb), provided the interactions are sufficiently strong. Furthermore, their 
analysis leads to the conclusion that persistent currents are unstable for I > 1 even when the concentration of the 
minority component is arbitrarily small. This latter conclusion seems at odds with what one might expect in the pure 
single-species limit (xb — > 0). 

In Sec. HVi we present the analysis of the GP energy functional in somewhat more detail than was provided 
by Smyrnakis et al. Q This analysis essentially confirms all of their analytical results, however, we find that the 
information regarding the behaviour of eo(L) in the vicinity of L = Na^ is not sufficient to establish whether or not 
persistent currents are actually stable. In fact, a more global analysis of e$(L) shows that persistent currents can 
exist when I > 1. Our work also clarifies how the single-species results are recovered in the xb — > limit. 

II. BLOCH'S CRITERION FOR PERSISTENT CURRENTS IN A TWO-SPECIES GAS 

In this section we extend Bloch's analysis to a two-species system consisting of Na particles of type A and Nb 
particles of type B. The masses of the particles are Ma and Mb- In addition, we assume an idealized one-dimensional 
ring geometry. The Hamiltonian H for this system is taken to be 



where the angular momentum operator of the i-th particle about the centre of the ring is 



The index i denotes an A- type particle for 1 < i < Na and a B-type particle for Na + 1 < i < Na + Nb = N. The 
subscripts on the interaction potential u,j allow for the interactions between the particles to be species-dependent. 
For the pair-wise interactions assumed, the total angular momentum 




(4) 




(6) 



commutes with the Hamiltonian. The stationary states ^(0i, •••,#jv) °f the Hamiltonian can thus be chosen to be 
simultaneous eigenstates of the total angular momentum. 
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A suitable basis of states can be constructed from the following product states for noninteracting particles: 

H9 1 ,...,9 N ) = 4> mi (e 1 )4> m2 (e 2 )---4> m!f ^N)- (7) 

Here rrii is an integer and 

p imO 

<j> m (6) = -=. (8) 
v Zir 



The wave function in Eq. ([7]) is an eigenfunction of L with eigenvalue L = h'^ li rrii. It can be written in different 
ways. One possibility is 



$(6>i,...,6W) = (2iry N/2 exp{iNl@) 



exp 



N 



(9) 



where 

A' 



1 N 

Z = -5> 4 (10) 



is the angular momentum per particle in units of h and 

1 N 

Q =N^ (ID 

i=l 

is the mean angular coordinate. The above wave function is identical in form to that of a single species system. 
By construction, the first exponential in Eq. @ is an eigenfunction of L with eigenvalue L = NIH. The second 
exponential is a function of the coordinate differences 0i — 9j and as such is a zero total angular momentum wave 
function. 

Properly symmetrized functions are obtained from Eq. © with the application of the symmetrization operator 

S = S A S B (12) 

where 

A ' Pa 

i '-mZ p '- (14) 

The operator Pa permutes the coordinates of the A particles, whereas Pb does the same for B particles. Applying 
the symmetrization operator to the wave function $(#i, On) yields 

^ {mi} {6 u ...,6 N ) = e^{iNlQ)x{ mi }{9u-,0 N ), (15) 

where X{m ; } is a normalized function of the coordinate differences 0i — 9j. The functions in Eq. (|15[) provide a basis 
of properly symmetrized Af-particle states, with N = Na + Nb- 

The stationary state solutions of -ffvl/ = E^> with angular momentum L will be denoted ^La(0i, On), where a 
indicates the rest of the quantum numbers. These states can be expanded in terms of the basis functions Eq. (|15|) as 



*Ia(9l,.-.^) = ^'CWW)*^}^!,---,^) 

{mi} 

= exp[iNie]xLa(6i,...,e N ), (16) 

where the prime on the summation implies the restriction XL =1 rrii — Nl. It is clear from the way Xha{9\, 9n) is 
defined that it is a function of the relative angular coordinates 9i — 0j . Substituting Eq. (|16|) into the Schrodinger 
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equation for $i a , we find that Xha satisfies the equation 

HlXLcc = e a (L)xLa, (17) 



where 



NR 2 \ M A M B I V ' 



and 



(19) 



- m Tp f t\ f h 2 l 2 N A h 2 l 2 N B \ L 2 (N A N B 

e a (L) - E a (L) - ^^-^ + J - E a (L) - [m^ + m^ 

Since L\La = 0, Hl in Eq. (TTTf can be expressed equivalently as 

HL=H+ wk(ik-i^) iA - (20) 

We observe that this Hamiltonian is in general L-dependent which has important consequences for e a (L). 

Eq. (|T7|l must be solved with appropriate boundary conditions. Since the wave function ^ La is required to be 
single- valued with respect to each of the angular variables, it satisfies 

9 La (...,0 i + 2*,...) = 9 La (...,0 i ,...). (21) 

Eq. (|2"Tj) then implies 

XLa (-, l + 2tt, ...) = cxp [-^2 7 r/] X L a (-, &i, •■•) (22) 

for i — 1, N. From this we see that the boundary conditions are periodic as a function of L with period Nh. With 
the basis functions written in the form given in Eq. ([9]), they are the same boundary conditions that apply in the 
single-species case. In the Nb = limit, e a (L) = E a (L) — L 2 /2MtR 2 = e a (L) with Mt = NMa- In addition, 
the Hamiltonian Hl reduces to H since La — > L and Lxlu = 0. As a result, XL'a with L' = L + Nh satisfies the 
same Schrodinger equation and boundary conditions as XLa ■ This implies that the eigenvalue spectrum for these two 
functions is identical. As concluded by Bloch [1], the eigenvalues e a (L) for the single-component system are then 
periodic functions of L with period Nh. In particular, the ground state energy is given by 

ML) = ^^ 2+ e (L). (23) 

The same considerations apply to the two-species situation for the special case Ma = Mb since Hl in Eq. (|2"0)) also 
reduces to H in this limit and Eq. (|2"3")l is still valid. In both of the above situations, the periodicity of the eigenvalues 
e a (L) means physically that the "internal" excitations can be the same for distinct macroscopic flows whose angular 
momenta differ by some multiple of Nh. However, eo(L) is no longer periodic when Ma ^ Mb, since the Hamiltonian 
Hl in Eq. (jTT)) depends explicitly on L. 

An alternative analysis is provided by writing the wave function in (|16j) as 



*ta(«i, On) = exp(iATO cm ) X La(0i, 9 N ), (24) 
where cm is the 'centre-of-mass' angular coordinate defined as 

1 N 



'cm 

Mr 

i—i 



and 



XLa(0i, ...,9 N ) = exp 



I N 



-i 

M T 



XL a (0i,...,e N ). (26) 
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Here and in the following, Mj is equal to Ma for i < Na and Mg for i > Na', Mt = NaMa+NbMb is the total mass. 
We observe that the exponential in Eq. (IM1) is still an eigenfunction of L with eigenvalue L and that Xha is a function 
of coordinate differences and therefore a zero-angular momentum function. Eq. (f24| amounts to a separation of the 
centre-of-mass motion from the internal degrees of freedom. Indeed, substitution of Eq. into the Schrodinger 
equation for yields 



HxLa = e a (L)xLa, (27) 

where 

e ot {L)=E a {L)-^- 2 . (28) 

Eqs. I|27p and (f2"8")l suggest that \La and e a (L) can be viewed, respectively, as the "internal" wave function and 
"internal" excitation energy. The boundary conditions imposed on XLa(9i, ■■■,0n) can be derived from Eq. (|24[) and 
are given by 

*£«(• • • , 9i + 2tt, • • • ) = exp ^27T-^^ X za(- • ■ , 6 U ■ ■ ■ ), (29) 

where v = Nl. When Ma = Mb — M, these boundary conditions revert to those of the single-species case where 

Xlc{- ■ ■ ,&i + 2tt, • • • ) = exp (-»27rl) xt«(- • • , ^i, • ■ • )• (30) 

This, together with Eq. (|27|) implies that X£+JVfi,a = XLq and e Q (L + NH) = e a {L). In fact, in this case e a (L) and 
Xha coincide with e n (L) and XLa, respectively. 

When Ma =/= Mb, e a {L) is not in general a periodic function of L. However, it can be if the boundary conditions 
in Eq. (|29p remain unaltered when v is augmented by some number N (i.e., L — > L + iVft) such that 

-mV =p > (31) 

and 

NM B 



Mt 



9, (32) 



where p and q are both integers. This implies that Ma/Mb must be equal to the rational number p/q. The lowest 
possible value of TV is obtained when p and q have no common divisor and is then given by 

N=pN A + qN B . (33) 

With this choice of N, e a (L) is a periodic function of L with periodicity Nh. In this situation, it is possible to 
impart a definite angular momentum to the two-species system without altering its "internal" state. For two different 
atomic species, the mass ratio Ma/Mb is never strictly a rational number and thus e a (L) cannot be strictly periodic. 
However, if 

Ma/Mb ~ P /q + 6 (34) 

where \S\ << p/q, one would expect e a (L), by continuity, to be quasi-periodic with a periodicity close to (NAP+Nsq)h. 
For example, a mixture of 85 Rb (A) and 39 K (B) has a mass ratio 

Ma/Mb ^ 2 + 0.07, (35) 

in which case the quasi-periodicity of e a (L) would be (2Na + Ns)fi. 

In the rest of this section we discuss the close connection between Bloch's argument on persistent currents and 
Landau's criterion for superfluidity. Our analysis mainly concerns the single-species and equal-mass two-species 
systems, where there is strict periodicity for e a (L). However, it also applies to the two-species system with unequal 
masses, insofar as it is a good approximation to regard &o{L) as quasi-periodic. According to Bloch, persistent currents 
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can occur at the angular momenta L n — nNh, for integral n, if Eq(L) has a local minimum at L = L n . We thus 
examine the behaviour of E (L) in the neighbourhood of L n . From Eqs. ([T]) and © one has 

, (L n + AL) 2 , . 
E {L n +AL) = 2MtR 2 + e °( £ » + A£ ) 

n n AL + E (AL) 7 (36) 



2M T R 2 



where 0„ = L n / (M^R 2 ) is the angular velocity of the centre of mass of the system at L n . This expression for the 
energy is analogous to the expression obtained via a Galilean transformation for a homogeneous system in which an 
excitation is produced in the rest frame of the superfluid [3]. To make this correspondence evident, we define the 
velocity v n = L n /MxR and write the energy in Eq. (151)1) as 

E {L n + AL) = \M T v 2 n + ( ^ ] v n + E (AL). (37) 



2 1 " v R , 

The first term on the right hand side is identified as the kinetic energy of the superfluid moving with velocity v n . 
Likewise, the last term is identified as the energy of a stationary superfluid containing an excitation with "momentum" 
AL/ R. It should be noted, however, that the analogy is not complete since for a homogeneous system the superfluid 
velocity v n can take arbitrary values whereas for the ring geometry the angular momentum is restricted to the discrete 
values L n . 

With this correspondence in mind, we take Eq(AL) to be the energy of the system with a single quasi-particle 
excitation with angular momentum AL = hm and energy e(m), i.e., 

E o (AL)=E o {0) + e(m). (38) 

We thus have 

E {L n + AL) = E (L n ) + e(m) + fl n hm. (39) 

The stability of the state with energy Eo(L n ) is then assured if the excitations lead to an increase in energy. In other 
words, the system will sustain persistent currents at L n for an arbitrary excitation of the system if 

e(m) + Ml n m > (40) 

for all m. Since e(— m) = e(m), the left hand side has a minimum for negative values of m and we thus require 

Vn<(pr) . (41) 

We have thus shown that Bloch's argument for persistent currents in the one-dimensional ring geometry naturally 
leads to the more familiar Landau criterion for superfluidity. If e(m) has a positive curvature as a function of m, 
which precludes a roton-like minimum, the inequality in Eq. (|41[) can be replaced by 



n n < (42) 

It is clear from this expression that the inequality must eventually fail when n exceeds some critical value n cr . 



III. BOGOLIUBOV EXCITATIONS, DYNAMIC STABILITY AND PERSISTENT CURRENTS AT 
INTEGER VALUES OF ANGULAR MOMENTUM PER PARTICLE 

The Landau criterion derived in the previous section focuses attention on the elementary excitations of the system. 
In this section, we obtain these excitations for a two-species gas in a one-dimensional ring geometry in the Bogoliubov 
approximation. We then apply Eq. (1421) to discuss persistent currents at integer values of angular momentum per 
particle for an equal-mass two-species system. 

In the following, we assume that the particles interact via contact interactions with strengths U ss > , where s, s' — A, B 
specify the species. Using the single-particle basis in Eq. the Hamiltonian in Eq. @ can be written in the second- 
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quantized form 



k = EE 



E E 



4~ 



i s,m u 's',Ti-m u 's',m'"'S,n-m' , 



(43) 



where to is the angular momentum quantum number and e s — h 2 m 2 /2M S R 2 . Assuming both species to be Bose- 
condensed in the to = state, the corresponding Bogoliubov Hamiltonian can be written as 



\ e s + 9ss)0, s . m a s,m + 2S ssa s,m a s,-m ' ~^9ssO'S,rnQ'S,-'m 



1 



t ~t 



1 



(44) 



where g ss , = U ss > y / N s N s >/2ir. 

The diagonalization of a Hamiltonian similar to Eq. (|44l) for a three-dimensional system was carried out in [f| . Here 
we present a different method of determining the Bogoliubov quasiparticle operators. This is done in three steps. 
First, we perform a Bogoliubov transformation for each of the species treated individually. The transformation is 
defined by 



with 



where 



*s,m a s,m 



a s,mt J s,m u s,ml J . 



(0) St 
s,m< s, — m 



<-s,— m — a s,ml J s,—m u s,mt J s,rm 



,(°) 



>2 1 /e 



egg - ^: 



(45) 



(46) 



E s = v/ e 2 + 2e s3;j 



(47) 



_E S is the Bogoliubov excitation energy for independent components. Substituting Eq. (|45j) into Eq. (|44[) and dropping 
all constant terms, we obtain 



^Bog = E E E Amkrn + £ £ 



(48) 



where 5 = t A ^B I E a Eb9ab- The second term in this Hamiltonian describes the coupling between the Bogoliubov 
excitations defined for each of the species. It is convenient to write the Hamiltonian (again to within a constant) in 
the matrix form 



#Bog = E 



(49) 



m>0 



where 



* m = A<m (3 A _ m fr B ,m /3s,_ m ) T 



and 



M = 



(E A g g 

E A g g 

9 9 Eb 

V g g E B 



(50) 



(51) 



To complete the diagonalization process we introduce the following transformations 

/3 S , -m = ^3^71^+,— m ~ Vs~jm@+,m ^s.m/^-.-m ~ ^s,m(^-,Tm (52) 

where the amplitudes are chosen to be real. The Hamiltonian is reduced to the diagonalized form 
if the amplitudes satisfy the matrix equation 

a z A4\v± = w±w± (54) 

with the normalization condition 

w±<7 z w± = 1. (55) 
Here, w± = {u^} n — U Bm — Vb^J t and the matrix cr z is defined as 



10 
0-100 
10 
-1 



(56) 



It should be noted that Eqs. (pyfj) and (|55|) guarantee that the Bose commutation relations of the new operators p+ m 

and /3-, m are preserved. 

The Bogoliubov excitation energies E± are determined by the characteristic equation 

det(cr z X - El) = 0, (57) 

which yields 

(E 2 - E\) (E 2 - E 2 B ) - 4e A e B g 2 AB = 0. (58) 
This quadratic equation in E 2 has the two roots [6j, l7j 



El = \(E\+E B )±\ y/(E A + E%f + 4 (4e A e B g 2 AB - E\E%). (59) 

The dispersion of these modes is 'phonon-like' for small m (E± cx |m|) and 'particle-like' for large m (E± cx m 2 ). The 
upper branch has the higher sound speed and evolves continuously into h 2 m 2 /2M < R 2 where -M< signifies the smaller 
of the two masses. 

The Bogoliubov excitations of the two-component system are dynamically stable provided E 2 > 0. Since only E 2 _ 
can become negative, the criterion for dynamic stability is 

E\E% - 4e A e B g AB > 0. (60) 

In view of Eqs. (f57| and (|58|) . this is equivalent to the condition 

dct(a z M) = det(M) > 0, (61) 

since det(cr z ) = 1. Using the definition of E 2 in Eq. (14T1) and defining 
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Eq. (|60|) becomes 

1 



iaa + 2 m J y BB + 2 m ) >lAB - ^ 

For repulsive interactions, this inequality is satisfied for all m if it is satisfied for m = 1. This limiting case gives the 
condition 

IAA + ~) (iBB + ^ > (64) 

A criterion of this form was obtained in |4( for Ma = Ms but is also seen to be valid for Ma ^ Mb with the definition 
of j ss i given in Eq. (|d^|) . 

To complete our discussion of the Bogoliubov excitations we present the results for the Bogoliubov amplitudes. It 
is straightforward to show that Eqs. and (|55j) lead to 

[u s , m ) ±E±E s {2El-E\-E%) 1 ' 

(~(±)\ 2 = (E±-E s )HEl-Ef) 

[Vs > m > 4E ± E S (2E 2 ± ~E 2 A ~ E%) ■ 1 ' 

where s denotes the species complementary to s. Finally, the relation of the original creation and annihilation 
operators to the Bogoliubov quasiparticle operators is defined via 

Q-s.-m — u i^~mft+,-m — v i^nfl+,m + u s,rnft-,-m — v i lT ny-,m- (67) 

These amplitudes can be obtained from Eq. (|46]l and Eqs. (|65l) and (l66t with the result 



(±)} 2 = (E ± +e s ) 2 (E 2 ± -E 2 ) 
4E ± e s (2E 2 ± -E 2 A -E 2 B ) 



(^±- £s ) 2 (^-g|) 



= ,r: ^ v "' 2 =gk - (69) 



It can be shown that these expressions are equivalent to those given in Ref. § in the one-dimensional limit. The 
amplitudes can be used to evaluate the mode density fluctuations Sni^n(6) of each species. We find that the A and 
B density fluctuations are in-phase for the (+) mode and out-of-phase for the (-) mode. 

We now make use of these results in Eq. PS|) to investigate the possibility of persistent currents at the angular 
momenta L n = nNH for the equal-mass system. The lower of the two branches in Eq. (|59|) is the branch relevant to 
determining the stability of the current. For Ma = Mb = M, the energy of this branch reads 



E-(m) = 2 ^ R2 ]jm A + m 2 (^/ A A + Ibb - \[{lAA - Ibb) 2 + ^Iab^J ■ 



(70) 



According to Eq. (|42j) . the stability of persistent currents at L n requires 



nh h / f ~~ 7, 

Jib 2 < 2MB 2 V 1 + 1AA + 1BB ~ V (7AA ~ 1BB) + AlAB ' (71) 

This inequality is satisfied if the following two inequalities 

Iaa g I [ lBB 2 J > 7j4s ' ^ 

Iaa + Ibb > 4n 2 -l, (73) 
are simultaneously satisfied. In the limit jab = 0, we have two independent components and we observe that the 
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inequalities are satisfied if 7 m i n = mia.(TAA, Tbb) satisfies 

Tmin > i(2n + l)(2n-l). (74) 

For n = 1 this gives the critical interaction strength 7 cr — 3/2 which is the value quoted in Ref. 

For the two-species system with equal masses, the inequalities in Eqs. (|T2[) and ([73|) can usually be satisfied for 
suitable choices of the interaction parameters, implying the possible stability of persistent currents at any L n . The 
only exception occurs when 

Iaaibb = Tab, ( 75 ) 

or equivalently 

UaaUbb = U 2 AB . (76) 



In this case, the coefficient of m in Eq. f|T0[) vanishes and the lower branch has a free particle dispersion which 
destabilizes persistent currents for any value of n. This conclusion was arrived at earlier by Smyrnakis et al. |4| for the 
special case V aa = Ubb = Uab] we see here how it follows from the Landau criterion for the more general relation 
in Eq. (|76|) . However, this does not preclude the possibility of persistent currents at non- integral values of angular 
momentum per particle. In the next section we reconsider the problem from the point of view of mean-field theory, 
following closely the work of Smyrnakis et al. Q 



IV. PERSISTENT CURRENTS AT NON-INTEGER ANGULAR MOMENTUM PER PARTICLE: 

MEAN-FIELD THEORY 



The analysis in this section is based on the mean-field Gross-Pitaevksii energy functional for the two-component 
system in the ring geometry: 



1* * (ass* 



dtpj 



d6 

271 



N B h 2 

2M B R 2 



dO 



2 U aaN a 



Here the condensate wave functions %p A and ip B are normalized as 



f d6\^ A \ 4 + lu BB N B [ d0\il> B \ 4 + U A BN A N B f d0\ijA\ 2 \i>B\ 2 - (77) 
Jo * Jo Jo 



d0|lte(0)| s 



(78) 



As discussed in the previous section, Bloch's argument allows for persistent currents at integral values of I = L/Nh 
when Ma = Mb = M except when Eq. (|76|) is true. Here, following Smyrnakis et al. [H, we consider the special case 
Uaa = U BB = Uab = U. In units of the energy Nh 2 /(2MR 2 ), Eq. §77} becomes 



f de(x A 


dipA 


2 


dips 


dB 


+ x B 


dO 





+ x 2 a ttt d6\ijj A \ 4 + x 2 b ttt d9\^ B \ 4 + 2x a x b ttt d6\iP A \ 2 \ijj B \ 2 , 



271 



(79) 



where xa = Na/N, x b = N B /N are the relative fractions of the two species in the system and 7 = NMR 2 U/irH 2 is 
a dimensionless interaction parameter. For definiteness, we take Na > N B . 

The objective is to minimize the energy functional in Eq. (|79p with the constraint that that the average value of the 
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total angular momentum has a fixed value L = INK. This is achieved by expanding the condensate wave functions as 

MO) =$>„^ m (0) (80) 

m 

^ B {0) = J2 d ^(0), (81) 

i n 

where the basis functions <j>m(0) are given in Eq. (J5J). The normalization of the wave functions requires 

]T| Cm | 2 = l, ^\d m \ 2 = l. (82) 

m tci 

Such a superposition implies that the wave functions are in general nonuniform around the ring. In addition, the 
expansion coefficients c rn and d m must satisfy the angular momentum constraint 

I = x A l A + x B h = xa ^m|c m | 2 + x B ^m|<i m | 2 . (83) 

m m 

I A (Ib) represents the average angular momentum in units of h of an A (_B)-species particle. The minimization of the 
energy with respect to the expansion coefficients in Eqs. (|80|) and (|8Tj) was first considered by Smyrnakis et al. Q. It 
will be clear from the following that much of our analysis closely follows theirs. However, we have expanded on their 
discussion in order to obtain a number of results that are not given explicitly in their paper. 
Substituting the wave functions in Eqs. (|80|) and (|8lT) into Eq. (|79|) . we obtain 

/-27T 4 

E {1) ^ x A ^m 2 \c m {l)\ 2 + x B ^m 2 \d m {l)\ 2 + x 2 ^ j d6^c m (l)(j> m (9)\ 

r^n . 4 />2tt . 2 

+ x|7T7 / d9\^2d m (l)(f) m (6) / ^^(i)^^) d m (l)<j) m (9) 

Jo m Jo m m 

= / 2 + e (/). 

According to Bloch's argument, eo(l) should exhibit the periodicity eo(l + n) = eo(l) where n is an integer. This 
periodicity is ensured if the expansion coefficients satisfy the periodicity conditions 

c m +n{l + n) = c m (l), d m+n (l + n) = d m (l). (84) 

The fact that Eq(1) must remain unchanged when the wave functions with angular momenta — l a are used to 

evaluate the energy functional leads to the relations 

c m (-l)=c*_ m (l), d m {-l) = d*_ m {l). (85) 

These two conditions are the mean-field counterparts of Eqs. ([2|) and ([3]). 

The function eo(l) is the central quantity determining the possibility of persistent currents and its detailed evaluation 
is taken up next. To begin, we consider wave functions ipA and %ps containing only two components, that is, 

■ipA = c (/)o + ci<t>i (86) 
ip B = d (fio + difa. (87) 



The coefficients c m and d m are normalized according to Eq. (|82j) and the angular momentum constraint becomes 

x A \ci\ 2 + x B \di\ 2 =1. (88) 

Expressing the complex coefficients in the form 

Cm = \c m \e ta ™ (89) 
d m = \d m \e i ^, (90) 
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the GP energy becomes 

Eo(l) = I + ~ +7 04M 2 |ci| 2 + x||d | 2 |di| 2 + 2s it ajfl|co||c 1 ||do||di| cos X ) (91) 
where \ — at® — a\ — /3q + j3\. The choice of x which minimizes Eq(1) is it and we then have 

Eo{l) =1 + ^+l{xA\c \\c 1 \-x B \d \\d 1 \) 2 . (92) 
The lowest possible value of this energy is [1] 

E (l) = l + j/2, (93) 

which occurs for 

za|co||ci| = a;.B|<io||di|. (94) 
This relation, together with the normalization and angular momentum constraints, yields the coefficients 

2 (x A -l)(l-l) 2 l(x B - I) , Q _v 

a;A(l — 2/) £yi(l — 2/) 

2 (^g ~ 0(1 - , ,2 1{XA-1) (Qf ,s 

x B (l-2/) x b (1-2/) 

These quantities are positive provided I is in the range < I < xg ot x A < I < 1. Assuming the validity of Eq. (|9"3")l 
for I in these ranges, we see that Eq(1) does not have a local minimum at I = l. Thus, persistent currents are not 
possible at I — 1, and by virtue of the periodicity of eo(l), at all integral values of I. These conclusions are consistent 
with our earlier discussion based on the Landau criterion; the validity of Eq. ([76| implies the existence of particle-like 
excitations and the absence of persistent currents at integral values of I. 

Although Eq. (|93[) was obtained for the simplest possible variational wave function, it in fact is exact when I is 
restricted to the above ranges To show this, we consider normalized wave functions of the form 

$a(0) = iPa{0) + S^ A , $ B (9) = ip B {0) + Sip B , (97) 

where ipA and ipB are defined by Eqs. (|86l) and (1871) with the coefficients given in Eqs. (|95|) - (|96)) . If the deviations are 
expressed in the form 

Sip A = X! Sc ^<p7n, SlpB = ^ Sd ^<Prn, (98) 

m m 

the angular momentum constraint in Eq. (|83|) leads to 

x A ci(Sci + 5c*) + XB<ii{8di + 8d*) — —xa J]] m|(5c T „ | 2 — xb /J m\5d m \ 2 . (99) 

rn m 

We next observe that the density n$(6) = Na\iPa\ 2 + Nb\iPb\ 2 is in fact uniform, that is, 710(6*) = N/(2ir). Using 
these results, the energy is found to be given by 

i-2-k 

E$A,1>b] =^o(0 + ^5Z(m 2 -™)|'5c m | 2 + a;B ^(m 2 ~m)| ( 5d„ l | 2 + ^ / d9\Sn(6)\ 2 , (100) 

where 8n(8) — n(6) — no. We thus see that E[i/ja, ips] > Eq(1), implying that the state defined by Eqs. (|86p and ([57| 
is indeed the ground state of the system for the assumed ranges of the angular momentum. It should be noted that 
this result depends crucially on the assumption of equal interaction parameters between all components. The weaker 
condition in Eq. (|76|) still precludes the possibility of persistent currents at integral values of I, but the energy does 
not have the simple form shown in Eq. (|93[) . 

We next analyze the energy for xb < I < xa- In particular we consider the situation when I is close to xa, that is 
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I - x A 



-e, where e is a small positive quantity. For I — xa we see from Eqs. (|95|) - (|96)) that 



\c \ 2 = 0, 
\d \ 2 = l, 



\ci\ 2 = l 
Mi| 2 = 0. 



(101) 
(102) 



As e increases from zero, we therefore expect deviations from these limiting values and additional components in the 
expansion of the tpA and i^b wave functions. To be specific, we consider the three-component wave functions 



ljj A = C (f>o + C\4>\ + C 2 (f)2 

tpB = d-i<j>-i + d cj) + dx(j)i. 



(103) 
(104) 



We anticipate that |co| 2 , | C2 1 2 , M-i| 2 and \d\\ 2 are all of order e. With this assumption, the energy to first order in e 
is found to be 



E (l) = x A (| Cl | 2 + 4|c 2 | 2 ) + x B (M_i| 2 + Mil 2 ) + 1 - 

+ a; 2 4 7 (|co| 2 + |c 2 | 2 + 2|c ||c 2 |cosxi) +x| 7 (|d_ 1 | 2 + |rf 1 | 2 + 2|rf_ 1 ||di|cosx 2 ) 



(105) 



2xaxbJ 



+ I cq 1 1 — 1 1 cos \3 + \co\\di\ cos(x3 - X2) + |c 2 ||d_i| cos(x3 - Xi) + MMil cos(x 3 - xi - X2) 



where we have defined the phase angles xi = a o — 2«i + ct 2 , \2 = P-i — 2/?o + Pi and X3 = a o — a i — + A)- This 
energy is an extremum with respect to the phase angles if they are all or 7r. If we choose them arbitrarily to be 0, 
we obtain 



E (l) ~ x A (|ci| 2 + 4|c 2 | 2 ) + x B (IcLxl 2 + K| 2 ) + |+7 MM + M) + a; B (|d_i| + |d x 



(106) 



which must now be minimized with respect to the coefficients |co|, |c 2 |, and |di| subject to the angular momemtum 
constraint 



I = x A - e = x A (| Cl | 2 + 2|c 2 | 2 ) + x B (l^l 2 - |d_x| 2 ) =x A (l- \c \ 2 + |c 2 | 2 ) + x B (|di| 2 - M-i| 2 ) 



(107) 



If this minimization in the end leads to coefficients that are negative, the phases have to be adjusted accordingly to 
yield coefficients with positive values. As we shall see, this will indeed be necessary. 

Using Eq. (|107l) to eliminate |ci| from Eq. (|106l) . and introducing a Lagrange multiplier A to account for the angular 
momentum constraint, the functional to be minimized is 



^(|c |, |CaU \ d -il Mil) = 2za|c 2 | 2 + 2x B M_ 1 | 2 + 7 x A (|c | + |c 2 |) + x B (M-i| + Mil) 



+ A 



X A (1 



CO 



C2 



X B (-M-ll 



Mil 



where the variations of the coefficients are now unconstrained. This variation leads to the results 



A + 2' 



-1, 



d-i 



co 



A-2' 



where the Lagrange multiplier A is the solution of the cubic equation [i[ 

/(A) = A(A 2 - 4) - 2 7 A + ^(xa - x B ) 



0. 



(108) 



(109) 



(110) 



The roots of this equation are to be determined for 7 > and < (xa — x B ) < L 

In Fig. [l] /(A) is plotted for (xa ~ xb) = 0, 0.5 and 1 and for 7 = 2; Figj2] is a similar plot for 7 = 8. For 
(xa - xb) = 0, /(A) = A(A 2 - 4 - 27), which has the roots A = and A = ±y/4 + 27. For (xa ~ xb) = 1, 
/(A) = (A - 2)[A(A + 2) - 27], which has the roots A = 2 and A = -1 ± y/1 + 27. The latter two values arc the 
Lagrange multipliers in the single-species limit as obtained from the minimization of Eq. f| 108[) for xb — 0. Since the 
term 4rf(xA — xb) in /(A) simply shifts the curves in Figs. Q] and [2] vertically, it is clear that there are always three 
real roots for the physical range of (x A — Xb) values. For any positive value of 7, one root is always less than —2, a 
second lies in the range < A < 2 (more precisely in the range < A < 2(xa — xb)) and a third in the range A > 2. 
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FIG. 1: Plot of the cubic /(A) vs. A. The curves from bottom to top correspond to xa — xb — 0, 0.5 and 1.0. The 

interaction parameter is 7 = 2. 
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FIG. 2: As for Fig. [I] but for an interaction parameter of 7 = 



Substituting the coefficients given in Eq. (|109j) into Eq. (|107j) we find 

e(A 2 -4) 2 



A[x A {\ + 1)(A - 2) 2 + .t b (A - 1)(A + 2) s 



(111) 



It is clear from this expression that the A < —2 root makes |co| 2 negative. This root is therefore physically inadmissible 
and only the positive A roots are relevant. Eq. (|lllj) together with Eq. (|109[) can be used in Eq. (|106p to evaluate the 
energy. One finds the remarkably simple result 



E o(l) - - = x A ~ £ + \e = x A + (I - x A ){l - A). 



(112) 



We now see that the smaller of the two positive A roots gives the lowest possible energy. This thus identifies the root 
in the range < A < 2 as the one that is physically relevant For A in this range we observe that the ratios in 
Eq. (|109[) are negative, indicating that the phases in Eq. (|105[) were chosen incorrectly. The proper phases are xi = ^ 
X2 = and X3 = 71". 

The criterion for the existence of persistent currents at I — xa used in Ref. [J] is that the slope of Eq(1) in Eq. (|1 12[) 
at I — x^ is negative, i.e., A > 1. The critical condition is thus A = 1, which from Eq. (jllOl) gives the critical 
interaction strength Q 



A(x A - x B ) 



2{Ax A - 3) ' 



(113) 



In the xa = 1 limit this reduces to 7 cr = 3/2 which is the value obtained at I = 1 for the single-species system. To 
obtain the critical coupling at / = x A + n — 1, where n = 1, 2,.., we write Eq(1) = I 2 + eo(l) and use the fact that eo(/) 
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FIG. 3: The angular momentum change carried by each of the wave function components relative to the total 
angular momentum change of —e as a function of xb- red (AIa,o), black (AZ^), green (AIb,-i), blue (AIb,i)- The 

interaction parameter is 7 = 2. 



is periodic. The slope at / = [xa + n — 1) is thus found to be 

dE (l) 



dl 



= 2n-l-\. (114) 

l=(xA+n — l)~ 



If the root in the range < A < 2(xa — Xb) is used, the slope cannot be zero for any n > 1. This is the basis of the 
claim made in Ref. [J] that persistent currents are not possible for I > 1; seemingly, an arbitrarily small amount of the 
minority component B has a profound effect on the possibility of persistent currents. For the single-species case, the 
energy is given by Eq. (11 12p with xa = 1, but the appropriate value of A is A = —1 + y/1 + 27, which is not bounded 
as a function of 7. Using this value in Eq. (|114j) . one finds that persistent currents are possible for all n in this case, 
with a critical interaction strength of j C r.n = (2n + l)(2n — l)/2. This is the result found earlier (Eq. ([74]) ) using the 
Landau criterion. This comparison indicates an inconsistency. On the one hand, Eq. (|114[) does allow for persistent 
currents for I > 1 in the single-species limit if the appropriate value of A is used. However, the two-species analysis 
requires that the root in the range < A < 2 be used, which precludes the possibility of persistent currents for I > 1 
for any nonzero value of xb- Since the energy functional in Eq. (|79[) reduces to the single-species case when xb = 0, 
it would appear that taking the xb — > limit of the two-species analysis is problematic. 

In order to explain this discrepancy it is useful to examine the behaviour of the coefficients in Eqs. (|109[) and (|111|) 
in the xa — > 1 limit in more detail. These coefficients are determined by the root A that lies in the range < A < 2. 
If 7 < 4, the limiting value of this root for xa — > 1 is A = — 1 + + 2j. This is the A value for the single-species 
case. Thus for this range of 7, one recovers the single-species values for all the coefficients. However, for 7 > 4, the 
root in the range < A < 2 has the limiting value of 2 which is less than the A = — 1 + y/1 + 2-f root. The limiting 
values of the coefficients do not correspond to the single-species values in this case. 

The distinction between 7 < 4 and 7 > 4 is revealed more clearly by plotting the coefficients in these two cases as 
a function of xb ■ We observe that the angular momenta carried by each of the species is given by 

Ia = x A (|ci| 2 + 2|c 2 | 2 ) =xa + x a (M 2 - |c | 2 ) (115) 

and 



l B =x B {-\d- 1 \ 2 + \d^ 2 ) . (116) 

The change in angular momentum as I is reduced from xa is associated with the transfer of weight from one angular 
momentum component to another. For example, for the A species, the transfer takes place from the m = 1 state to 
the ui — or vn — 2 states, with respectively, a decrease or increase in angular momentum. For the B species, the 
transfer takes place from the m = state to the m = — 1 and m — 1 states. Of interest is the relative magnitude 
of the angular momentum change Al s>m that is associated with each angular momentum component. We therefore 
define the ratios Al s<m /(— e) where for example, AIa,o/{^£) = (— xa\cq\ 2 )/(— e). These ratios represent the fraction 
of the angular momentum change —e — I — xa attributable to each of the angular momentum components. In Fig. [3] 
we plot these ratios as a function of xb for 7 = 2; Fig. @] gives similar plots for 7 = 8. For 7 = 2, we see that species 
B carries a relatively small contribution of the angular momentum change. This contribution vanishes in the xb — > 
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FIG. 4: As in Fig. 1 but for 7 = 8. 



limit and the situation reverts to that of the single species which, as discussed above, is generally the case for 7 < 4. 
The situation for 7 > 4, however, is quite different. Fig. |4] for 7 = 8 shows that the angular momentum change is 
carried entirely by the m = — 1 component of the B species in the xb — > limit. The reason for this surprising 
result is that the relevant A root approaches 2 for xb — > when 7 > 4. Eq. then gives | Co | 2 — (2 — \) 2 e/(4xb) 

and from Eq. (|109j) we find |gLi| 2 — s/xb for xb — > 0, i.e. lg = — e in this limit. The divergence of |<i_i| 2 as 
xb — > is indicating that the result can only be valid for a decreasingly smaller range of e since the normalization 
1 = |rf-i| 2 + \d \ 2 + \di\ 2 must be preserved. In other words, the energy Eq(1), as given by Eq. (|112|) . is meaningful 
in an interval of / of decreasing size as xb — > 0. 

The above results call into question any conclusion regarding the possibility of persistent currents at higher angular 
momenta when x a approaches 1. In this limit, a more global perspective regarding the behaviour of the energy 
as a function of I in the interval xb < I < xa is required. We now give a general argument for the possibility 
of persistent currents at I > 1 based on the assumption of continuity of the GP energy as a function oi xb- To 
exhibit this dependence we write Eq(1,xb) and consider this function in the limit of small xb- In particular, we 
have E (n,x B ) = E^(n) + 5%(xb) and E (n — Al,x B ) = E^(n — Al) + 5 2 (x B ) where Eq(1) = E (l,x B = 0) is the 
energy of the single-species system. The assumption of continuity implies that 5\{xb) and ^(^b) approach as 
xb — > 0. We then have E (n — AI,xb) — E (n,XB) — Eq(u — Al) — Eq(ti) + ^(^s) — $i{xb)- By choosing 7 > 7 cr ,n, 
Eq(ti — Al) — Eq(u) will have some fixed positive value. Thus, we can say that Eq(ji — AI,xb) — Eo(n,XB) > for 
xb sufficiently small. Since 8Eq(1, Xb) /dl\i =n - < 0, we conclude that Eq(1,xb) must have a local minimum between 
I = n — Al and / = n. This argument can be used for any n and shows that persistent currents must be stable in the 
vicinity of I — n if xb is sufficiently small and 7 is sufficiently large. 

Although it is difficult to evaluate Eq(1,xb) for arbitrary I, the above general argument can be illustrated quanti- 
tatively by evaluating the energy at I — 1/2. To do so, it is sufficient to assume four-component wave functions of the 
form 



*I>A = C-\<j>-\ + C Q (p Q + C\$\ + C 2 (f> 2 

i) B = d-i<f>-i + d Q 4> + di<j>i + d 2 4>2- 



(117) 
(118) 



Substituting these wave functions into Eq. (IT9"|) . we have 



E (l = 1/2) = x A (|c_i| 2 + | Cl | 2 + 4|c 2 | 2 ) + x B (M-i| 2 + Mil 2 + 4|d 2 | 2 ) 





+ x A Ti-fj dQ\c-i4>-i + c Q if)Q + cx<f>i + c 2 (f>2\ +x B ir~f d0\d-i(j>-i + d cf)o + di<fri + d 2 4>2 \ 
Jo Jo 

r 2ir 



(119) 



The periodicity and reflection properties imply cq (|) = c\ (i) and c_i (|) = c\ (i), with analogous relations for the 
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d m amplitudes. These relations reduce the number of variational parameters by half. We have in particular 



|co| = |ci| = x, 
ot x = -a , 

\do\ = \dx\ = u, 
Xi = -Xo, 



\c-i\ = |c 2 | = y 
a_i = -a 2 
\d-i\ = \d 2 \ = v 
X-i = -X2- 



Using these definitions, the normalization constraints reduce to 



x 2 + y 2 



1 2 

2' 



v 2 = -. 
2 



Furthermore, the angular momentum of each species is given by 

l A = x A {-\c-i\ 2 + | Cl | 2 + 2 1 c 2 1 2 ) = x A {x 2 + y 2 ) 
l B = x B {~\d-i\ 2 + Mi| 2 + 2|d 2 | 2 ) = x B {u 2 + v 2 ). 

We thus sec that normalization ensures that the total angular momentum has the required value of 1/2. 
Using these results, the expression for the energy becomes 

£ (l/2) = \ + \l + Ax AV 2 + 4x B v 2 



(120) 



(121) 



(122) 
(123) 



x* + y 4 + 8x 2 y 2 + Ax 3 y cos j3 



x 2 B l 



u 4 + v 4 



8u 2 v 2 



Au v cos £ 



+ X A X B 1 



Sxyuv {cos(6» - /3 + £) + cos(26> - /3 + £)} + Axyu 2 cos{9 - 0) + Ax 2 uv cos{9 + £) 
+ 2x 2 u 2 cos 9 + 2y 2 v 2 cos(36» - 2/3 + 2£)1 , 



(124) 



where we have defined the phase angles /3 = 3«o + <*-2, £, = %Xo + X2 and — 2(ao — Xo)- We see that the energy 
depends on these three independent phases and the two amplitudes x and u. It clearly reduces to the single-species 
result in the x B — > limit. 

For xa = 1, the energy is minimized for j3 = n and a value of x which is close to l/y/2.. We do not expect this 
conclusion to change when x A is close to, but not exactly equal to 1. For these values of xa, the term in Eq. (|124[) 
proportional to x 2 B is small and can be neglected. Setting /? = 7r, the energy is approximately 

S (1/2) ~ ~ + \p + Ax A y 2 + Ax B v 2 + x 2 Al x 4 + y 4 + 8x 2 y 2 - Ax 3 y 



+ X A XBl 



- 8xyuv {cos(# + £) + cos(2# + £)} - Axyu cos 9 + Ax uv cos(9 + £) 
+ 2x 2 u 2 cos 6> + 2y 2 v 2 cos(36< + 2£)1 , 



(125) 



From this we see that the phases 9 and £ only appear in the last term proportional to xb- It is clear that Eq is 
stationary with respect to these phases when they take the values and it. To explore the various possibilities, we 
define the function 



f(x, u, £, 0) = - 8xyuv[cos{9 + £) + cos(26> + £)] - Axyu 2 cos 9 

+ 4a: 2 m; cos(6> + £) + 2x 2 u 2 cos 9 + 2y 2 v 2 cos(36» + 2£) , 



(126) 



which is the quantity multiplying x^^b! hi Eq. (|125|) . This function is tabulated in Table U for various values of £ 
and 9. From this table it is clear that £ = 0, 9 = it will give a lower energy than £ = ir, 9 = tt. For £ = 0, (9 = we 
have 



f(x, u, 0, 0) — f(x, —u, 0, 0) = 8xuv(x — Ay). 



(127) 



Since xa is close to 1, Eq. (|125[) is minimized for a value of x close to l/y/2 which is much larger than y. This 
implies that any minima of the function f(x, u, 0, 0) will occur for negative values of u. But u must be positive (recall 
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f(x,u,£,9) 








-IGxyuv — <ixyu + Ax uv + 2x u + 2y v 





7T 


Axyu 2 — Ax 2 uy — 2x 2 u 2 — 2y 2 v 2 


TT 





+ I6xyuv — 4xyu 2 — 4x 2 uv + 2x 2 u 2 + 2y 2 v 2 


TT 


7T 


Axyu 2 + Ax 2 uy — 2x 2 u 2 — 2y 2 v 2 



TABLE I: The function f(x,u,£,9) defined in Eq. (|126[) tabulated for various values of £ and 9. 




FIG. 5: The energy at I = 3/2 vs 7 for xa — 0.95. The horizontal line is the value of Eq{1 + x A ) — 7/2. 



u = \do\), so this case must be rejected. Finally, for £ = 71", = 0, we have 

/(a;, w, 7r, 0) — /(x, — u, 7r, 0) = — 8xuv(x — Ay). 



(128) 



The same argument implies that minima of f(x,u, 7r, 0) must occur at positive u. We are thus left with the two 
possibilities £ = 0, 9 = ir and £ = 7r and 6> = 0. A comparison of the contour plots of Eo(x, u, ir, 0) and Eq(x, u, 0, it) 
shows that the latter is the one that provides the lowest energy. For xa — 0.95 and 7 = 2, Eq(x, u, 0, ir) is minimized 
for x m in ~ 0.697 and u m i n ~ 0.677. The value of x m i n found here is close to the value of 0.696 found for xa = 1. Not 
surprisingly, the |c m | 2 coefficients are close to the values obtained in the single-species limit. 

We will now use the value of Eq(1/2) to show that persistent currents are possible for I > 1. To be specific, we 
consider 1 = 1 + 1' with < /' < 1. Using the periodicity of eo(0; we have 



£?o(l + = 1 + 2/' +E (l'). 



(129) 



At I = x A , Eq. ([93j gives E (x A ) =x A + 7/2. We then find that E (l + x A ) -7/2 = 1 + 3x A = 3.85 for x A = 0.95. 
As explained earlier, this value is exact within the mean-field analysis. We next use Eq. fj 129[) to obtain 



E (3/2) = 2 + E (1/2) . 



(130) 



In Fig. [5] we show the behaviour of Eq(3/2) ~ 7/2 as a function of 7 for xa = 0.95. We see that -Eo(3/2) becomes 
larger than £?o(l-95) at a value of 7 ~ 15. This implies the existence of a local minimum in the range 1.5 < I < 1.95 
and hence the possibility of persistent currents. The value 7 ~ 15 is clearly an upper bound to 7 cr for this value of 
x A - 

The approximate behaviour of Eq(1) as a function of / can be obtained by generating approximations to e~o{l)- For 
< I < x B and x A < I < 1, e (Z) - 7/2 = Z(l - I). From Eq. (|112p we have e' \ l=x - = 1 - 2^ - A. The simplest 
approximation to e~o(l) in the range xb < I < a^A consistent with this information is 



e o 1} (0-7/2 = /(l-0 + A 



(x A - I) (I - x B ) 
xa - x B 



(131) 
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FIG. 6: The function eo(l) — 7/2 plotted vs I in different approximations. The blue curve is the function 1(1 — I); 
the red, black and green curves are and respectively. 7 = 16 and xa — 0.95. 



An improved approximation is a fit that reproduces the value of eo(l) at I = 1/2. It takes the form 



4 2) (i)-7/2 = Kl-0+ A 



(x A -l)(l- x B ) , {x A ~l) 2 (l - x B f 



XA ~ %B 



(x A - x B y 



(132) 



where fi = 16[eo(l/2) — 7/2 — 1/4 — X(xa — xb)/4], A third approximation ignores the information about the slope 
of eo(l) at I — xa but includes the value at I = 1/2. This approximation gives 



4 3) (o- 7/2 = i(i-i)+ v 



(x A - I) (I - x B ) 
(xa - xb) 2 



(133) 



where v — 4(eo(l/2) — 7/2 — 1/4). These various approximations are plotted in Fig. [6] for 7 = 16. We expect the 
correct variation of Eq(1) to be bounded by the &q (I) and (I) curves; for / — > 0.95, eo(l) should be closer to the 

f2l (3) (3) 

e (I) curve but for / — > 0.5 it should be closer to the (I) curve. We note that e (I) must give the correct 
behaviour in the xa — > 1 limit. 




FIG. 7: The energy E (l) 



7/2 vs I for 7 = 16 and xa — 0.95. The various curves correspond to the various 
approximations to &o(l) shown in Fig.[Sl 



These different approximations can be used to determine corresponding approximations to Eq(1), which is plotted 
in Fig. [7]m the range < I < 2 for 7 = 16. The red curve based on gg (I) does not show a local minimum at / = 1.95 

as predicted by considerations of the slope of E (l) at this point. On the other hand, the black curve based on Sq (I) 
which includes the information about eo(l/2) shows a local minimum below / = 1.95 and demonstrates that persistent 
currents should be possible for / between 1.5 and 2. Regarding persistent currents at I = 0.95, the critical interaction 
strength according to Eq. (|113j) for xa — 0.95 is j CI ~ 1.9. For this value of 7, eo^(l/2) is actually quite close to 
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to the true value eo(l/2) as determined by the four-component wave function analysis. Thus the prediction of the 
critical interaction strength based on the slope of E (l) remains quite accurate in this case. However, it is clear that 
the slope of Eo(l) calculated at I = (xa + n — 1) — , although correct, is not sufficient to provide a criterion for the 
existence of persistent currents for n > 1. We have also extended the analysis to slightly smaller values of xa and 
arrive at similar conclusions. However, increasingly larger values of 7 are then required to achieve a local minimum 
between / = 1.5 and 1 = 1 + xa- 

We finally mention the behaviour of Eq(1 = 1/2) when xa = 1/2. In this limit the expression for Eg (I) given in 
Eq. (|93l) is correct for all I and gives in particular Eq(1 = 1/2, x a = 1/2) = 1/2 + 7/2. This value is reproduced by 
Eq. (|124[) at xa = 1/2 irrespective of the phases (3 and £ since the minimum occurs for 6 = n and x — u — l/y/2, 
where all the j3 and £ dependent terms have no effect since y = v — 0. We note that the minimizing value of 9 is the 
same as in the xa ~> 1 limit and anticipate that this will remain true for intermediate values of xa between xa = 1/2 
and xa = 1- However, a more careful analysis of Eq. (|124[) would be required to confirm this and to determine the 
remaining variational parameters that minimize the GP energy. 



V. CONCLUSIONS 



In this paper we have extended to the two-species system Bloch's original argument regarding the possibility of 
persistent currents in the idealized one-dimensional ring geometry. Strict periodicity of the energy e a (L) defined in 
Eq. (JTJ is found to arise when the mass ratio Ma/Mb is a rational number. By making a connection to the Landau 
criterion for the special case Ma = Mb, we show that persistent currents are in general possible at the discrete set 
of total angular momenta L n = nNH, except when the interaction parameters satisfy the condition in Eq. (|76[1 . The 
underlying reason for this limitation is the existence of excitations with a particle-like dispersion. This conclusion is 
consistent with the predictions of a mean-field analysis based on the GP energy functional. A detailed analysis of 
the GP energy in the vicinity of I = xa, first carried out by Smyrnakis et al. [4j, indicates that persistent currents 
are possible at this angular momentum per particle if the interaction parameter exceeds the critical value given in 
Eq. (I113p . These authors go on to claim that persistent currents cannot arise for / > 1 in the two-species system. 
However, a more detailed analysis of the global behaviour of the GP energy demonstrates that this conclusion is 
not valid. Quite generally, the properties of the two-species system evolve continuously to those of the single-species 
system as the concentration of the minority component is reduced. It would of course be of interest to verify these 
theoretical predictions experimentally. The recent experimental realization of toroidal Bose-Einstein condensates [1, [s| 
would suggest that experiments on two-species systems may soon be feasible. 
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